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PREFACE 

This  report  was  prepared  by  Paul  V.  Sellmann,  Geologist,  of  the  Northern  En- 
gineering Research  Branch,  Experimental  Engineering  Division  of  the  U.S.  Army 
Cold  Regions  Research  and  Engineering  Laboratory  (USA  CRREL),  Dr.  Jerry 
Brown,  Research  Soil  Scientist,  Earth  Sciences  Branch,  Research  Division,  USA 
CRREL,  R.l.  Lewellen,  Arctic  Geologist  (P.O.  Box  1068,  Littleton,  Colorado 
80120)  and  Dr.  Harlan  L.  McKim,  Soil  Scientist,  and  Carolyn  J.  Merry,  Geologist, 
Earth  Sciences  Branch,  Research  Division,  USA  CRREL.  The  study  was  performed 
under  DA  Project  4A161 102B52E,  Research  in  Military  Engineering  and  Construc- 
tion, Task  02,  Research  in  Snow,  Ice  and  Frozen  Soil,  Work  Unit  014,  Environ- 
mental Factors  Influencing  Expedient  and  Permanent  Construction  of  Logistic 
Facilities  and  Army  Operations  in  Cold  Environments. 

Technical  reviewers  of  this  report  were  Dr.  Robert  F.  Black  of  the  University 
of  Connecticut  and  Richard  Haugen  of  USA  CRREL. 

The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or 
promotional  purposes.  Citation  of  trade  names  does  not  constitute  an  official 
endorsement  or  approval  of  such  commercial  products. 
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THE  CLASSIFICATION  AND 
GEOMORPHIC  IMPLICATIONS  OF  THAW  LAKES 
ON  THE  ARCTIC  COASTAL  PLAIN,  ALASKA 

P.V.  Sellmann,  J.  Brown,  R.I.  Lewellen, 

H.  McKim  and  C.  Merry 

INTRODUCTION 

A major  portion  of  the  Alaskan  Arctic  Coastal  Plain  is  underlain  by  ice-rich  marine 
sediments,  the  product  of  several  marine  transgressions  and  regressions.  Numerous  thaw 
lake  basins  of  variable  morphology  and  distribution  have  developed  on  the  perennially 
frozen  sediments  (permafrost)  of  this  low-lying  plain.  Most  notable  are  the  large  elliptical 
and  oriented  lakes  found  on  the  northernmost  part  of  the  Coastal  Plain.  These  lakes  have 
long  been  considered  a unique  feature  of  this  environment  and  their  origin  and  geometry 
have  captured  the  interest  of  many  investigators  (Black  1969a  and  b.  Black  and  Barksdale 
1949,  Britton  1967,  Carson  and  Hussey  1962,  Livingstone  et  al.  1958).  The  following 
features  make  these  lakes  an  important  key  to  understanding  the  history  and  current 
processes  controlling  the  geomorphology  of  this  and  other  circumpolar  permafrost- 
dominated  terrains. 

1.  The  thaw  lake  cycle.  The  thaw  lake  cycle  (Britton  1967),  consisting  of  repetitive 
stages  of  lake  formation  and  drainage,  represents  the  primary  geomorphic  process  that 
modifies  this  land  surface.  Thaw  lakes  are  largely  unstable  with  active  erosion  at  their 
basin  margins.  As  a result,  they  often  coalesce  and  drain.  Therefore,  in  some  areas  where 
this  process  is  active,  these  lakes  seldom  exceed  several  thousand  years  in  age  (Brown 
1965,  Carson  1968).  Recently  a relatively  restricted  mechanism  of  thaw  lake  formation  was 
recognized  whereby  sediment  laden  ice  rafts  initiated  thawing  of  the  permafrost  and 
formation  of  lake  basins  (Ugolini  1975). 

2.  Topographic  conditions.  Topographic  conditions  control  lake  size  and  shape.  The 
lower  the  regional  slope  and  the  local  relief  the  greater  potential  for  lakes  of  large  size. 

3.  Erosion  and  local  relief.  Lake  shape  and  orientation  are  linked  to  two  principal 
controls.  The  m^jor  mechanism  in  areas  of  low  relief  and  ice-rich  sediment  is  preferential 
lake  erosion  induced  by  waves  and  currents  (Carson  and  Hussey  1962).  Another  control  is 
local  relief  resulting  from  geological  structure  (fault  lines,  basins,  lineaments)  or  the  topog- 
raphy associated  with  dunes,  periglacial  landscapes,  and  fluvial  patterns  (Fischer  and 
Lathram  1973,  Short  and  Wright  1974,  Gatto  and  Anderson  1975). 

4.  Ground  ice  distribution.  Thaw  lakes  influence  ground  ice  distribution  in  a region. 
Through  the  erosional  process  and  subsequent  freezeback  of  lake  sediments,  thaw  lakes 
control  and  tend  to  reduce  the  amount  and  type  of  near-surface  ground  ice  and  con- 
sequently the  surface  relef  (Hussey  and  Michelson  1966). 

5.  Basin  geometry.  In  general,  lake  depth  is  fairly  uniform  and  is  usually  less  than  2 m. 
These  lakes  commonly  freeze  to  the  bottom.  The  lakes  with  some  persistent  ice  in  July 
generally  correlate  with  the  deeper  lakes  (>  2 m)  that  usually  do  not  freeze  to  the  bottom 
(Brewer  1958).  Lake  depth  is  related  to  ground  ice  volume,  local  relief,  and  basin  age. 
Many  lakes  have  a similar  basin  geometry  with  a shallow  marginal  shelf.  The  deepest 
central  part  of  the  basin  is  usually  elongate  (Black  and  Barksdale  1949). 

Lake  depth  determines:  I)  the  feasibility  of  using  the  lakes  for  winter  water  supply,  2) 
the  biological  resources  associated  with  them,  and  3)  the  influence  of  the  lakes  on  the 
thickness  of  the  thaw  zone  beneath  and  adjacent  to  the  basins.  Only  the  lakes  that  do  not 
freeze  to  the  bottom  are  useful  for  water  supplies  and  contain  overwintering  fish  popula- 
tions. The  depth  to  permafrost  beneath  a lake  is  closely  related  to  the  depth  of  water. 
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APPROACH  AND  METHODS 


Earlier  studies  of  lake  distribution  and  morphology  have  been  based  on  topographic 
maps,  conventional  aerial  photography,  field  studies,  isotopic  dating  and  comparison  with 
other  coastal  areas.  With  the  advent  of  high  quality,  summer  Earth  Resources  Technology 
Satellite  (ERTS-I)  imagery,*  it  was  possible  to  view  the  regional  patterns  of  Alaskan 
Coastal  Plain  lakes  at  a scale  which  provides  simultaneous  examination  of  lake 
morphology,  distribution  and  sequential  changes  in  ice  cover.  Analyses  of  sequential  ice 
cover  changes  proved  to  serve  as  a valuable  indicator  of  water  depth.  In  addition,  land- 
water  area  data  previously  available  on  the  protracted  (BLM)  survey  maps**  have  been 
computer  digitized  and  the  results  computed  and  summarized  at  the  township  level.  This 
information,  in  combination  with  ERTS-1  imagery  and  our  familiarity  with  the  field  condi- 
tions and  genetic  processes,  permits  us  to  offer  additional  information  on  the  occurrence 
and  distribution  of  arctic  thaw  lakes,  and  their  geomorphic  and  historical  significance.  Our 
approach  was  both  a qualitative  and  quantitative  assessment  of  lake  morphology  and 
related  depth  parameters. 

A lake  classification  was  developed  which  categorized  the  lakes  based  on  basin 
morphology  and  distribution  patterns.  An  ERTS-1  mosaic  at  a scale  of  l:50Q,000  (Fig.  lb) 
was  compiled  using  1973  summer  imagery  (Multispectral  Scanner  System;  MSS  band  7, , 
0.8-1. 1 /am)  and  provided  the  working  base  for  this  study.  The  parameters  considered  in 
the  classification  include  size,  axis  orientation,  degree  of  basin  development,  and  lake 
density.  The  initial  grouping  was  based  on  visual  contrasts  and  patterns  seen  on  the 
mosaic  (Table  I).  Measurements,  from  standard  USGS  topographic  maps  (1:250,000  scale) 
of  over  1500  lakes  for  dimensions  and  basin  orientations,  provided  a quantitative  basis  for 
the  classification  (Fig.  2)  and  confirmed  the  more  subjective  grouping  established  by  visual 
comparison  from  the  ERTS-1  mosaic.  All  lakes  with  major  axes  0.50  km  or  larger  were 
measured.  Measurements  of  the  length  of  nu(jor  and  minor  axes  and  the  bearing  of  mtyor 
axes  were  obtained  from  the  USGS  topographic  maps  instead  of  the  ERTS-1  imagery  for 
several  reasons:  1)  the  maps  provided  a uniform  basis  for  determining  the  shape  or  limits 
of  the  lake  basins  which  the  ERTS-1  scale  and  sequential  images  did  not  provide;  2)  the 
resolution  of  the  maps  is  better  than  the  ERTS-1  imagery  for  measurements  since  they 
were  constructed  from  high  resolution  aerial  photographs  of  larger  scales  than  ERTS-1;  3) 
the  maps  provided  coverage  of  areas  that  were  partially  or  completely  cloud  covered  on 
the  ERTS-I  summer  mosaic.  The  tabulated  and  plotted  data  (Fig.  2)  were  developed  from 
selected  portions  of  each  lake  unit  and  covered  approximately  15%  of  the  entire  Coastal 
Plain.  The  number  of  observations  made  to  establish  these  values  in  each  lake  unit  is  also 
indicated  in  Figure  2. 


RESULTS 

Lake  morphology 

Size . Although  lake  size  varies  greatly  across  the  Alaskan  Arctic  Coastal  Plain,  a sub- 
division based  on  size  was  feasible  (Table  I)  using  three  classes:  large,  intermediate,  and 
small.  Lake  size  is  controlled  or  influenced  by  ground  ice  volume,  local  and  regional  relief, 
and  age.  The  amount  of  ground  ice  is  related  to  the  depositional  and  freezing  histories  and 
the  texture  of  the  material  in  any  specific  area,  finer  textures  yielding  large  quantities  of 
ground  ice.  The  largest  thaw  lakes  would  be  expected  on  flat-lying  surfaces  of  fine-grained 
sediment  which  contain  large  amounts  of  ground  ice  and  where  lake  expansion  is  uncon- 
strained by  local  relief. 

* The  imagery  used  for  this  study  was  obtained  from  NASA,  as  part  of  CRRLI.  contract  S-702S3-AG  KRTS-I 
Data  User  Investigation  of  Arctic  and  Subarctic  Environmental  Analysis  Utilizing  ERTS-I  Imagery:  D.M. 
Anderson,  principal  investigator  (l)b  329).  The  liRTS  satellite  has  subsequently  been  renamed  LANDSAT. 
On  file  in  the  Anchorage  office  of  the  Bureau  of  Land  Management. 
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ERTS-l  uncontrolled  mosaic  (a  composite  of  1973  band  7,  summer  imagery). 
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c.  Contour  of  percentage  land  cover  by  water 


Figure  I . Maps  of  northern  Alaska  showing  most  of  the  Coastal  Plain,  the  northern  foothills  and  portions  of  the 
northeastern  Brooks  Range. 
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Table  I.  Classification  and  major  features  of  the  lakes  of  the  Alaskan  Arctic  Coastal  Plain. 

Development  of  Orientation  of  Lake  density 

elonXaRe  Q-xis elongate  axis  (7,  water  coverage) 
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Figure  2.  Statistical  characteristics  of  lake  units  1-4. 


/ 


The  largest  lakes  occur  on  the  outer  Coastal  Plain  margin  in  the  L,  unit  and  exceed  6 km 
in  length.  This  is  an  area  of  low  regional  and  minimum  local  relief  as  compared  to  the  other 
parts  of  the  Coastal  Plain.  This  area  also  consists  of  fine-grained  marine  sediments  which, 
based  on  experience  at  Barrow  (Brown  and  Sellmann  1973),  have  a high  ground  ice 
volume.  Since  the  largest  lakes  occur  on  the  lowest  relief  portion  of  the  Coastal  Plain  (in 
the  Teshekpuk  Lake  area)  and  the  smallest  lakes  occur  in  the  L3  unit  (which  contains 
surfaces  with  greater  regional  relief)  an  attempt  was  made  to  test  the  hypothesis  that  there 
is  a maximum  lake  size  for  any  given  regional  slope.  A simple  step  approach  superimposed 
on  the  regional  slope  was  used.  This  approach  appears  valid  since  the  preferred  orienta- 
tion, or  maximum  length,  of  the  lakes  is  parallel  to  the  slope  of  most  surfaces.  The  only 
assumptions  required  for  this  approach  were  estimates  of  lake  scarp  heights  commonly, 
found  on  the  Coastal  Plain.  The  estimated  scarp  heights  used  were  1,  3 and  6 m.  These 
values  limit  lake  length  for  any  given  regional  slope.  The  limit,  being  the  contact  of  lake 
ends  as  the  lakes  increase  in  length,  is  simply  expressed  as 

Lmax  = h cot  a (1) 

where  Lmax  is  maximum  lake  length  in  meters,  h is  lake  scarp  height  in  meters,  and  cot  a is 
the  cotangent  of  the  regional  slope  angle. 

Regional  slopes  were  measured  and  maximum  lake  lengths  on  these  slopes  were  plotted 
on  both  logarithmic  and  arithmetic  graphs.  The  calculated  relationships  between  slope  and 
size  are  shown  on  the  arithmetic  plot  for  the  various  scarp  heights  mentioned  in  Figure  3. 


Q Regional  Slope  (degree*) 


Figure  3.  Lake  length  plotted  against  regional 
slope,  with  calculated  curves  superimposed  for  var- 
ious scarp  heights. 
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The  measured  regional  slope-lake  length  data  taken  from  standard  USGS  topographic 
maps  (1:250,000  scale)  cluster  in  the  range  of  the  scarp  heights  selected.  The  arithmetic 
plot  implies  that  regional  slope  has  a strong  control  on  the  potential  development  of  large 
lakes.  The  lakes  longer  than  5 km  can  only  occur  in  extremely  low  relief  areas  of  less  than 
0.1°,  and  as  the  slope  decreases  from  this  value,  lake  length  could  theoretically  approach 
infinity.  This  plot  also  indicates  that  at  low  regional  slope  angles  (s  0.7°),  lakes  longer  than 
0.5  km  should  be  rare.  It  therefore  appears  that  the  large  numbers  or  clusters  of  lakes  that 
are  apparent  on  1:250,000  scale  maps  and  ERTS-1  imagery  are  restricted  to  surfaces  with 
very  low  regional  relief.  In  areas  where  there  is  little  local  relief,  regional  relief  can  have  a 
major  control  on  lake  size. 

Elongation 

Another  common  feature  of  the  Coastal  Plain  lakes  is  their  development  of  an  elongate 
axis.  This  is  made  more  noticeable  since  most  of  these  elongated  axes  also  have  the  same 
preferred  orientation.  The  degree  of  development  of  elongate  axes  varies  usually  on  a 
regional  scale.  Based  on  visual  comparison,  the  variations  were  divided  into  four  groups: 
well-developed,  moderately-developed,  poorly-developed,  and  undeveloped  elongate 
axes.  The  measured  values  of  both  major  and  minor  axes  and  the  resulting  length/width 
ratios  provide  a check  on  the  above  qualitative  subdivision  of  these  parameters.  The 
average  length/width  ratio  for  the  unit  with  the  most  elongate  or  “well  developed”  elonga- 
tion was  2.1  and  the  “moderate  to  poor”  unit  had  an  average  length/width  ratio  of  1.6 
(Table  I and  Figure  2).  The  potential  for  axis  elongation  appears  to  be  strongly  controlled 
by  material  type  and  local  relief. 

Orientation 

In  this  arctic  coastal  setting  orientation  is  closely  related  to  the  degree  of  axis  elonga- 
tion, since  the  wind-induced  erosional  processes  which  encourage  elongation  and  orienta- 
tion are  the  same.  It  has  been  demonstrated  by  use  of  aerial  photographs  taken  over  a 
20-year  period  that  present-day  erosion  is  occurring  at  the  lake  extremes  (Brown  1970, 
Lewellen  1972).  These  observations  confirm  Carson’s  observations  that  elongation  occurs 
at  the  end  of  lakes,  partially  due  to  erosion  produced  by  wind  generated  waves,  currents 
and  higher  water  temperatures.  The  degree  of  orientation  for  classification  purposes  was 
divided  into  four  groups:  strong,  preferred,  poor,  and  no  orientation. 

All  but  one  of  the  main  lake  units  demonstrate  strong  or  preferred  orientation.  The 
control  of  orientation  is  much  the  same  as  for  elongation.  The  lack  of  orientation  can 
generally  be  related  to  shallow  overburden  leading  to  subsurface  control  or  surface  topo- 
graphic features,  such  as  dunes.  Rose  diagrams  of  long  axis  orientation  show  the  patterns 
for  the  various  classes  (Fig.  2). 

Percentage  cover  (density) 

the  estimate  of  lake  coverage  was  initially  based  on  a visual  comparison,  subsequently 
quantified  and  contoured  based  on  percentages  of  water-land  derived  from  the  BLM 
protracted  survey  maps  (Fig.  lc).  These  values  were  derived  from  acreage  values  for 
water  and  land  at  the  section  level  for  the  area  north  of  the  Umiat  Meridian.  A computer 
program  was  employed  to  calculate  the  percentage  of  water  per  section  and  make  summa- 
tions compiled  at  the  township  level.  Since  we  were  only  interested  in  the  fresh  water  on 
land  surfaces,  townships  with  saline  lagoons,  bays,  and  other  coastal  water  features  were 


6 


rejected  in  the  computer  summations.  The  final  township  values  were  manually  trans- 
ferred to  the  township  base  map  and  contoured.  Examination  of  Figure  lc  indicates  that 
the  valufes  for  total  area  covered  by  lakes  are  lower  than  might  be  anticipated  from  a visual 
inspection  of  maps,  aerial  observations  and  photographs.  For  example,  Hussey  and 
Michelson  (1966)  estimated  that  50-75%  of  the  Coastal  Plain  was  covered  either  by  lakes  or 
by  marshes  that  occupy  low  areas  or  former  lake  basins.  However,  this  study  showed  that 
only  very  few  areas  (not  exceeding  three  townships  in  size)  have  lakes  covering  more  than 
40%  of  their  surface,  and  these  areas  amount  to  only  about  1%  of  the  total  area  contoured 
in  Figure  lc  (although  marshes  and  drained  lake  basins  were  not  included  in  these  compu- 
tations). Approximately  only  22%  of  the  region  studied  has  lake  coverage  greater  than 
20%.  This  area  lies  on  the  outer  Coastal  Plain  and  is  characterized  by  surfaces  which  have 
not  been  modified  by  fluvial  processes.  The  wide,  alluvial  surfaces  have  a lake  coverage  of 
commonly  less  than  20%,  which  is  less  than  adjacent  Coastal  Plain  areas. 

The  larger  percentages  of  lake  coverage  are  found  on  areas  of  lowest  relief  on  the 
youngest  seaward  part  of  the  Coastal  Plain.  As  the  regional  slope  increases,  the  percent- 
ages of  water  coverage  decrease.  This  is  particularly  noticeable  on  the  high-relief,  inner 
margins  of  the  Coastal  Plain  where  values  are  below  5%.  Some  small  areas  of  larger 
percentages  are  localized  on  small  plateaus  or  low  relief  plains  isolated  in  areas  of  higher 
relief  (Fig.  lc). 

In  the  higher  relief,  western  section  of  the  Coastal  Plain  this  general  trend  is  again  borne 
out  where  low  percentages  of  water  cover  are  found,  as  compared  with  the  higher  percent- 
ages found  on  the  low  relief,  central  part  of  the  plain. 

In  some  areas  slope  may  only  be  a secondary  factor  controlling  the  amount  of  land 
covered  by  lakes,  with  such  factors  as  the  thickness  of  sediment  over  bedrock  being  more 
critical.  Although  many  areas  on  the  Coastal  Plain  have  very  small  contrasts  between 
slopes,  the  average  slope  is  still  large  enough  to  cause  major  differences  in  lake 
morphology.  Slope  is  assumed  to  be  the  primary  controlling  variable  in  differentiating 
units  L,  and  L3. 


LAKE  CLASSIFICATION 

Based  on  the  parameters  of  size,  development  and  orientation  of  elongate  axis,  and 
percentage  coverage  already  described,  the  lakes  were  grouped  and  mapped  into  six  lake 
units  and  three  nonlake  units  (Table  I).  Within  the  lake  units,  the  first  four  (L,-L<)  are  the 
most  common  and  well-defined.  The  remaining  two  (L5  & Lu)  are  more  general  and 
accommodate  the  few  departures  from  the  main  units.  Several  lake  units  contain  subordi- 
nate groups,  which  are  denoted  by  a letter  subscript,  and  represent  departure  from  charac- 
teristics within  the  main  unit. 

L,  unit 

The  Li  unit  is  characterized  by  large  lakes  which  cover  the  low-lying,  seaward  portion 
of  the  Coastal  Plain.  The  combination  of  low  regional  and  local  relief  and  the  properties  of 
the  sediments  are  thought  to  be  responsible  for  the  lake  basin  characteristics.  The  lakes 
are  most  distinctive  because  of  their  large  size,  elongated  shape,  and  dominant  long-axis 
orientation.  The  unique  appearance  of  the  lakes  in  this  unit,  in  contrast  to  those  in  the 
remainder  of  the  region,  suggested  to  early  workers  that  the  history  of  this  area  was 
different  than  that  of  adjacent  areas  (O’Sullivan  1961).  O’Sullivan  considered  this  unit  as 
the  youngest  part  of  the  Coastal  Plain,  which  must  have  been  formed  or  modified  by  a 
marine  transgression  i.i  late  Sangamon  or  mid-Wisconsin  time. 


Approximately  75%  of  the  lakes  in  this  unit  have  their  elongate  axes  oriented  between 
340°  and  360°.  This  is  also  the  only  unit  containing  lakes  that  significantly  exceed  6 km  in 
length,  and  this  unit  has  the  least  number  of  lakes  in  the  smaller  size  range  (1  km).  The 
lakes  have  the  largest  length/width  ratio,  supporting  the  visual  impression  that  these 
basins  are  more  elongated  than  lakes  in  other  units. 

A subordinate  group  LWa  was  established  for  a small  part  of  the  study  area  that  con- 
tained large  lakes  that  departed  from  the  main  unit.  Lakes  in  group  L,_a  have  poor  de- 
velopment of  elongate  axes  and  consequently  no  definite  orientation. 

L2  unit 

The  L2  unit  includes  lakes  that  are  large  to  intermediate  in  size.  None  of  these  lakes 
exceed  6 km  in  length,  and  the  number  of  smaller  lakes  in  the  1-km  size  range  is  greater 
than  in  L,.  The  basins  in  this  unit  are  more  elongated  than  those  in  all  of  the  other  units 
except  L,.  The  orientation  of  the  long  axes  of  the  basins  shows  the  same  general  strong 
trend  as  lakes  in  the  first  unit,  although  there  is  slightly  more  scatter  around  this  dominant 
orientation. 

Two  subordinate  groups  Lj-,  and  L2_t,  were  established  in  this  unit  to  cover  areas  in 
which  lakes  were  smaller,  less  dense,  or  less  oriented  than  in  the  mqjor  unit.  The  L2Hj 
lakes  were  primarily  found  on  the  higher  relief  inner  margin  of  the  Coastal  Plain,  and 
lacked  strong  orientation. 


L3  unit 

The  L3  unit  was  established  to  include  those  areas  covered  with  small  lakes  that  have 
strong-preferred  orientation.  The  histogram  of  major  axis  length  for  this  unit  shows  an 
increase  in  the  number  of  basins  in  the  1-km  size  range  and  no  lakes  with  axes  larger  than  3 
km.  The  rose  diagrams  of  axes  orientation  for  L3  lakes  show  the  same  general  trend  as  in 
the  previous  units. 

One  subordinate  group  Ls-a  was  recognized  for  small  areas  which  had  less  well- 
developed  orientation  and  lower  density  than  the  primary  unit.  This  unit  also  appeared  to 
contain  more  drained  lakes. 

L4  unit 

The  L4  unit  provides  a grouping  for  the  lake  basins  that  do  not  have  the  well-developed 
orientation  and  elongation  common  to  the  other  units.  These  lakes  appear  to  have  some 
control,  such  as  local  relief,  influencing  their  development.  In  some  cases  this  influence  is 
most  likely  related  to  surface  features  such  as  dunes  or  structural  control  in  areas  of  thin 
surface  sediment.  The  lengths  of  m^jor  axes  of  these  lakes  group  primarily  in  the  1-km  size 
range,  and  compared  to  the  other  lake  units,  this  unit  falls  between  L2  and  L3.  There  are 
also  few  large  lakes,  only  a very  small  percentage  of  which  exceed  6 km  in  length,  and 
compared  to  Li  and  L2,  this  unit  contains  a smaller  percentage  in  the  3 to  6 km  size  range. 
The  rose  diagrams  for  this  unit  show  the  pronounced  deviation  from  the  common  elonga- 
tion and  orientation  common  to  the  previous  three  units. 

An  additional  subunit  L«.a  was  established  to  accommodate  lakes  with  no  orientation 
and,  in  some  cases,  smaller  size  than  the  main  unit.  Although  not  shown  on  Figure  1 
because  of  size  limitations,  this  unit  is  the  common  one  on  the  eastern  section  of  the 
Coastal  Plain  in  the  vicinity  of  Barter  Island. 
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Ls  and  L„  units 

No  statistical  sampling  was  made  of  the  lakes  in  these  units.  These  lakes  often  lack 
preferred  orientation  and  elongation  of  their  basins  and  they  commonly  are  of  low  density. 
They  often  occur  on  the  inner  margins  of  the  Coastal  Plain  in  areas  of  high  relief.  Appar- 
ently the  sediments  and/or  the  relief  found  in  these  areas  have  an  overriding  control  on 
their  shape,  which  departs  from  the  characteristics  seen  in  other  parts  of  the  Coastal  Plain. 

Other  units 

The  remaining  units  primarily  contain  nonlake  surfaces  or  sections  of  the  Coastal  Plain 
on  which  lake  geometry  has  been  modified  by  more  dominating  processes,  such  as  fluvial, 
colluvial,  and  alluvial  activities. 

Alluvial  unit  (A).  The  alluvial  areas  are  usually  associated  with  current  river  systems, 
particularly  with  recent  flood  plains,  and  form  belts  or  strips  that  cut  across  the  lake  units. 
Several  of  the  A unit  areas  are  rather  extensive. 

Fluvial-lacustrine  unit  (F/L).  This  unit  includes  areas  of  low  relief  which  are  associated 
with  flood  plains  and  terraces  in  which  stream  channel  development  and  ground  condi- 
tions have  influenced  potential  lake  development.  It  is  most  common  on  the  lower  reaches 
of  the  more  major  rivers  and  streams  that  flow  across  the  flat  low-lying  portion  of  the 
Coastal  Plain.  Lakes  in  this  unit  commonly  have  a truncated  appearance  and  some  basins 
fill  old  meanders. 

Colluviallalluvial  unit  (C/A).  This  unit  was  developed  to  contain  what  appeared  to  be 
more  featureless  sloping  surfaces  along  streams  and  the  coastline  (areas  in  which  thaw 
lakes  are  not  common)  than  found  in  other  units.  This  unit  forms  bands  along  the  slopes  of 
more  major  streams,  and  the  transitions  between  land  surfaces  are  of  varying  elevations. 
The  slope  of  these  surfaces  may  be  as  low  as  2°  to  4°  based  on  observations  by  Hussey  and 
Michelson  (1966). 

Lake  basin  depths 

Depth  classes.  As  a result  of  the  examination  of  the  ERTS-1  imagery  for  development  of 
the  lake  classification,  additional  information  becomes  apparent,  lake  depth  being  the 
most  notable  and  significant.  Earlier  workers  had  indicated  that  lakes  with  a persistent 
mid-  to  late-summer  ice  cover  are  the  deepest  lakes  (Brewer  1958,  Black  and  Barksdale 
1949)  because  they  have  the  thickest  ice  mass  to  melt.  The  ERTS-1  imagery  enabled  us  to 
follow  progressive  melting  of  lake  ice  during  the  late  spring,  and  early-  and  mid-summer 
periods.  Using  all  ERTS-1  cycles  and  both  MSS  band  5 (0.6-0. 7 ^m)  and  band  7 (0.8-1. 1 
/am),  a fairly  detailed  inventory  and  classification  of  lake  depths  is  possible.  Three  depth 
categories  were  delineated:  1)  lakes  in  excess  of  2 m which  do  not  freeze  to  the  bottom,  2) 
more  shallow  lakes  of  intermediate  depth  (1-2  m)  that  freeze  to  the  bottom,  and  3)  very 
shallow  lakes  and  ponds  (<  1 m)  some  of  which  are  basins  that  are  seasonally  filled  at  the 
time  of  spring  runoff. 

The  approach  used  in  distinguishing  between  the  lakes  in  the  three  depth  classes  was 
based  on  June  and  July  1973  ERTS-1  imagery.  The  deepest  lakes  retain  their  ice  cover 
(white  on  ERTS-1  band  7)  for  the  longest  period.  During  June  and  early  July  an  open  water 
boat  is  formed  around  the  floating  ice  mass  which  corresponds  to  a shallow  shelf,  usually 
less  than  1 m in  depth,  which  surrounds  the  deeper  central  basin.  The  floating  ice  mass 
itself  corresponds  to  thicker  ice  formed  over  the  deeper  basin.  For  example,  a low  altitude 
aerial  oblique  view  of  a thaw  lake  in  early  July  shows  a large  ice  mass  remaining  over  the 
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deeper  central  part  of  the  basin  and  corresponding  in  shape  to  part  of  this  inner  basin  (Fig. 
4).  The  formation  of  the  shallow  moat  around  the  deeper  central  basin  can  be  seen  when 
comparing  Lake  A on  Figures  5a  and  b.  The  shallow  shelf  area  is  ice-free  early  in  July. 
Comparison  of  the  June  imagery  (Fig.  5a)  with  the  July  imagery  (Fig.  5b)  illustrates  the 
loss  in  ice  cover  during  the  elapsed  period. 


Figure  4.  A thaw  take  basin  containing  a shallow  marginal  shelf  and 
deep  central  basin  which  retains  a portion  of  its  ice  cover  in  the  deepest 
part  of  the  basin  (photograph  by  R.K.  Haugen,  USA  CRREL). 

The  second  depth  class  is  somewhat  less  obvious  than  the  first  group.  This  class  in- 
cludes lakes  that  have  frozen  to  the  bottom.  In  the,ERTS-l  MSS  band  7 imagery,  these 
lakes  appear  dark,  suggesting  no  ice  cover.  However,  this  class  of  lakes  appears  white  in 
MSS  band  5.  This  difference  is  due  to  freestanding  water  on  the  surface  of  ice  which  is 
frozen  to  the  bottom  and  the  greater  ability  for  water  penetration  with  MSS  band  5 
imagery.  A number  of  examples  of  this  can  be  seen  in  the  lakes  adjacent  to  Sinclair  Lake  in 
Figure  6.  A comparison  of  MSS  bands  5 and  7 imagery  of  Sinclair  Lake  also  suggests  that 
the  western  margin  of  this  lake  is  flooded  by  surface  runoff.  The  penetration  of  light 
through  water  for  MSS  band  5 imagery  has  been  reported  to  be  as  great  as  3 m (Polcyn  and 
Lyzenga  1973).  Brewer  (1958)  also  noted  that  the  ice  on  shallow  lakes  may  melt  from  the 
top  without  ever  floating,  and  that  this  ice  condition  can  frequently  be  observed  during 
flights  over  the  Coastal  Plain  in  late  June. 

The  third  depth  class  consists  of  the  shallow  lakes  that  apperr  dark  on  both  MSS  bands 
5 and  7 imagery  early  in  the  melt  season.  These  are  considered  very  shallow  lakes  (<  1 m) 
from  which  the  ice  cover  rapidly  melts.  They  also  may  be  only  seasonally  filled  lakes  that 
occupy  old  drained  lake  basins  at  time  of  runoff.  Some  of  these  seasonal  lakes  may 
completely  drain  or  the  water  may  evaporate  from  these  marshy  areas  by  mid-summer. 
Recently  Seilmann  et  al.  (1975)  appear  to  have  demonstrated  the  ability  of  side-looking 
airborne  radar  (SLAR)  to  differentiate  arctic  coastal  lakes  which  were  frozen  and  not 
frozen  to  the  bottom. 
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Ice  volume  and  basin  genesis 


The  depth  of  thaw  lakes  basins  appears  to  be  controlled  principally  by  the  ice  volume  in 
the  original  frozen  sediment  (Livingstone  et  al.  19S8,  Hussey  and  Michelson  1966).  Brewer 
(1958)  indicates  that  the  lakes  in  the  Barrow  area  fall  into  two  depth  classes:  0. 6-0.9  m and 
1 .8  to  2.7  m.  These  depth  ranges  can  be  considered  representative  of  the  most  outer  part  of 
the  Coastal  Plain. 

The  relationship  between  segregated  ice  volume  and  depth  was  developed  from  both  our 
studies  in  the  Barrow  area  and  from  those  of  O’Sullivan  (1966).  The  CRREL  data  are 
plotted  in  Figure  7 and  are  based  on  a composite  of  shallow  auger  holes  and  deeper  coring 
(Brown  1969,  Sellmann  and  Brown  1965,  1973).  The  noticeable  feature  of  this  curve  is  the 
large  amount  of  segregated  ice  in  the  upper  3 to  5 m.  Ice  volume  associated  with  ice 
wedges  is  not  included  in  this  plot.  For  purposes  of  potential  thaw  settlement  calculations 
(Fig.  8a  and  b)  estimates  of  wedge  ice  were  applied  to  the  upper  3 m of  ice  thickness,  based 
on  10  and  20%  additional  ice  due  to  wedges.  These  estimates  seem  reasonable,  based  upon 
observed  occurrences  of  wedges  associated  with  polygon  troughs  and  estimates  of  buried 
ice  (UNESCO  1970). 
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Figure  7.  Plot  of  ice  volume  vs  depth  for  the 
CRREL  auger  and  core  samples  obtained  in 
the  Barrow  area.  These  values  are  for  segre- 
gated ice  ( pore  and  lens)  and  do  not  include 
the  volume  ofic'  associated  with  ice  wedges, 
which  may  be  as  high  as  10  to  20%. 
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The  settlement  for  a given  amount  of  thaw  is  based  on  the  volume  of  ice  in  excess  of  the 
anticipated  volume  of  voids  formed  in  the  resulting  thawed  sediment.  The  excess  ice 
volume  was  directly  equated  to  potential  thaw  settlement  for  1-m  increments  and  plotted 
against  thaw  depth.  Two  porosity  values  corresponding  to  a silt  and  a sandy  silt  were  used 
for  the  thawed  sediment.  The  0.40  porosity  value  resulted  in  curves  that  flattened  at  a thaw 
depth  of  6 m,  since  at  the  7-m  depth  the  void  volume  exceeded  the  total  ice  volume, 
indicating  no  additional  potential  for  settlement.  The  curves,  based  on  a porosity  of  0.35, 
yield  a small  amount  of  excess  ice  at  depth  and  a potential  for  small  amounts  of  additional 
settlement.  The  real  situation  probably  falls  between  these  two  cases. 

The  data  indicate  a noticeable  decrease  in  ice  volume  with  depth  which  is  most  likely 
characteristic  for  most,  if  not  all,  of  the  Coastal  Plain  for  the  following  reasons:  1)  the 
limited  thickness  of  the  fine-grained  overburden,  which  has  the  potential  for  high-ice 
volumes,  overlying  low  ice  volume  bedrock  or  sediments,  and  2)  the  higher  volumes  of 
wedge  ice  and  lens  ice  which  are  more  commonly  found  near  the  surface.  Massive  ice  can 
be  found  at  greater  depths  under  special  conditions  particularly  under  rapidly  aggrading 
surfaces  (Mackay  1972). 

The  control  of  basin  depth  (depth  below  original  or  primary  surface  to  the  lake  bottom) 
on  the  data  in  Figure  8 becomes  more  obvious  when  it  is  appreciated  that  lakes  less  than  2 
m in  depth  freeze  to  the  bottom.  In  the  Barrow  area  these  shallow  ponds  and  lakes  have  a 
summer  thawed  zone  of  between  0.3  and  0.7  m which  refreezes  annually  (Brewer  1958, 
Coyne  and  Kelley  1974).  Lakes  deeper  than  2 m develop  a perennial  thaw  bulb  (Lachen- 
bruch  et  al.  1962).  Once  the  higher  ice  zone  has  been  thawed,  the  influence  of  the  progres- 
sively deeper  thaw  appears  to  be  small,  since  the  ice  volumes  with  depth  are  relatively  low 
and  thus  thaw  does  not  contribute  to  settlement. 

In  lakes  2 m deep,  bottom  temperatures  are  such  that  the  freezeback  potential  of  a 0.7-m 
annual  thaw  zone  is  marginal,  representing  the  threshold  at  which  a permanent  thaw  bulb 
will  start  to  form  (Brewer  1958). 

Depths  of  several  lakes  in  the  Barrow  area  were  plotted  on  thaw-settlement  curves  (Fig. 
8b).  Lake  depth  and  thaw  depth  is  also  presented  in  Figure  8c,  with  a being  distance  from 
the  original  surface  to  the  lake  surface,  b representing  maximum  lake  depth,  and  c indicat- 
ing depth  of  thaw  under  the  basin.  Therefore,  the  sum  of  a and  b equals  thaw  settlement, 
and  the  sum  of  a,  b and  c equals  total  thaw  depth.  Data  for  lakes  1,  4,  5 and  6 were 
obtained  from  observations  by  Brewer  (1958).  Although  limited,  actual  lake  data  appear  to 
agree  with  the  thaw-settlement  data;  no  thaw  depth  datac  were  available  for  lakes  5 and  6, 
indicated  by  dashed  lines  across  the  plot  (Fig.  8b). 

The  solid  line  x in  Figure  8a  indicates  that  shallow  lakes  have  no  potential  of  becoming 
deeper  than  1.4  m (a  + b)  in  terrain  containing  10%  by  volume  wedge  ice  and  a seasonal 
thaw  of  0.7  m.  Any  additional  thaw  settlement,  as  indicated  by  the  intersection  of  this  line 
and  the  thaw-settlement  curves,  would  require  additional  excess  ice  volume. 

A 20%  volume  of  wedge  ice  is  necessary  for  developing  basins  that  exceed  2 m depth. 
The  deeper  range  of  lakes  (4, 5 and  6)  in  the  Barrow  area  also  indicates  that  once  this  2-m 
range  is  exceeded  only  minimal  settlement  takes  place.  Lake  4,  for  example,  had  a 58-m 
section  with  temperatures  above  0°C  below  the  basin  and  only  3.3  m (o  + b)  of  settlement. 

In  their  investigation  of  tundra  relief  features  Hussey  and  Michelson  (1966)  discussed 
the  influence  of  thaw  consolidation  on  the  development  of  relief  on  the  Coastal  Plain  near 
Barrow.  In  their  approach,  they  considered  the  influence  of  thaw  lake  history  on  the 
various  surfaces  and  how  it  modifies  ground  ice  volume  distribution.  They  found  that  the 
oldest  surface  that  has  not  been  recently  modified  by  thaw  lakes  has  the  greatest  potential 
for  thaw  consolidation.  Their  grouping  of  the  land  surfaces  in  order  of  decreasing  ground 
ice  volume  is  as  follows:  initial  surface,  ancient  drained  lake  basins,  recent  drained  lake 
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Lake  Data 

(m)  Settlement  Thaw  Depth 


a.  The  solid  pair  of  curves  was  obtained  by  using  a 
wedge  ice  volume  of  10%  in  the  upper  3 m and  porosity 
values  for  the  resulting  thawed  sediment  of 0.35  and  0.40. 
The  other  pair  is  for  20%  wedge  ice  in  the  upper  3 m. 


1 a 0.2  (a+b)  (a+to+c) 

b 0.8  1.0  1.7 

c 0.7 

2 a 0.6 

b 0.9  1.4  2.1 

c 0.7 

3 a 0.2 

b 2.0  22  3.0 

c 0.8 

4 a 0.3 

b 2.8  3.1  61.1 

c 58.0 


5 a 0.6 

b 1.8  2.3 

c 

6 a 0.5 

b 2.7  3.2 

c 


c.  Sample  lake  data,  with  a being  the  distance  from  the  original  surface,  b repre- 
senting maximum  lake  depth,  and  c indicating  depth  of  thaw  under  the  basin. 


b.  Lake  data  (from  8a  and  8c)  plotted  for  comparative  purposes. 


Figure  8.  Potential  thaw  settlement  or  total  basin  depth  vs  total  thaw  depth  for  lake 
basins  using  ground  ice  volume  data  from  Barrow. 
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basins,  and  present  lake  areas.  The  ranges  of  settlement  for  these  surfaces,  based  on 
limited  subsurface  data  due  to  thaw  of  the  upper  6 m of  the  permafrost,  were  3.4  m,  1.2  m, 
0.7  m,  and  0.1  m respectively. 

Our  thaw  settlement  curves  provide  information  to  greater  depth,  since  they  are  based 
on  a large  number  of  drill  holes  from  all  the  above  surfaces  except  present  lake  basins. 

The  thaw  lake  settlement  process  for  the  Barrow  landscape  is  depicted  in  Figure  9,  an 
idealized  cross  section  between  Barrow  Village  and  Elson  Lagoon.  A series  of  lakes  has 
occupied  most  of  the  area.  The  highest  areas  in  the  Barrow  vicinity  are  considered  primary 
or  initial  land  surfaces  (Hussey  and  Michelson  1966).  However,  small  lakes  have  locally 
thawed  and  reworked  some  of  this  surface.  This  primary  surface  is  the  source  of  much  of 
the  ice  volume  data  in  Figure  7 to  depths  approaching  sea  level.  Immediately  to  the  east  of 
this  area,  a series  of  old  lakes  has  thermally  eroded  the  original  ice-rich  sediments  to  about 
5 m above  sea  level,  some  3 m lower  than  the  adjacent  primary  land  surface.  This  amount 
of  thaw  settlement  agrees  with  Hussey  and  Michelson’s  determinations  as  well  as  with  our 
own  thaw  settlement  curve.  Furthermore,  radiocarbon  dating  of  a buried  lacustrine  peat, 
indicating  an  age  of  12,160  ± 200  years  B.P.  (Lewellen  1972)  at  2.2  m above  sea  level, 
suggests  the  maximum  period  over  which  this  process  has  occurred. 

The  basin  shown  in  Figure  4 illustrates  the  characteristic  shallow  shelf  surrounding  a 
deeper  inner  basin  common  to  some  lakes  in  this  setting.  The  origin  of  these  shelves  is  not 
completely  understood.  Black  (1969b)  indicates  that  the  shelves  may  have  a complicated 
history  in  that  portions  of  them  may  be  depositional  and  other  parts  erosional.  It  is 
suggested  by  the  authors  that  the  shelves  of  some  lakes  may  be  closely  linked  to  the 
thermokarst  origin  of  these  lakes,  and  very  likely  to  a particular  stage  in  their  development 
cycle.  Adequate  ground  ice  must  be  available  in  the  permafrost  to  provide  local  settlement 
of  the  basin  in  excess  of  2 m upon  thaw.  As  indicated  earlier,  when  this  critical  lake  depth 
is  reached  in  the  central  part  of  the  basin,  a perennial  thaw  zone  will  develop  beneath  it 
(Fig.  10).  The  increasing  thaw  depth  can  result  in  additional  settlement  in  the  central  basin 
if  adequate  ground  ice  is  available.  On  shallower  shelves  only  a seasonal  thaw  zone  is 
maintained  with  no  additional  thaw  settlement  due  to  the  insulating  effect  of  accumulating 
sediments  and  shallow  water.  The  thaw  of  the  permafrost  underlying  the  central  basin  may 
eventually  continue  toward  the  lake  margin  until  most  of  the  basin  is  deepened  to  a depth 
greater  than  2 m.  This  would  represent  a more  mature  stage  in  lake  development  in  areas 
where  subsurface  conditions  did  not  limit  basin  development. 


b.  Control  Batin  (>2  meters  In  depth) 

Figure  10.  Idealized  cross  sections  demonstrating  how  the  interior 
basin  of  a shallow  lake  (a)  can  deepen  through  thaw  settlement  as 
the  lake  exceeds  the  critical  depth  of  2 m (b). 
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GEOLOGICAL  IMPLICATIONS 


The  ERTS-1  imagery  and  its  regional  view  of  the  arctic  lakes  established  not  only  the 
influence  of  local  and  regional  relief,  and  to  a lesser  degree  ground  conditions  on  lake 
geometry,  but  also  provided  further  evidence  on  the  historical  development  of  the  Coastal 
Plain.  In  developing  the  depositional  history  of  this  region,  past  workers  have  relied  in  part 
on  the  physical  appearance  of  the  land  surface,  mainly  exemplified  by  the  lake  distribution 
patterns  and  topographic  breaks.  The  positions  of  previously  recognized  strand  lines  are 
shown  in  Figure  1 1 . When  comparing  the  positions  of  these  strand  lines  it  can  be  seen  that 
those  suggested  by  Black  (1964)  and  O'Sullivan  (1961)  are  in  general  agreement  for  the  last 
major  transgression  during  mid-  to  late-Wisconsin  time.  This  limit  coincides  with  the 
surface  covered  by  large  lakes  in  the  L,  unit.  Some  variation  in  lake  types  exists  in  this 
region,  largely  caused  by  disruption  of  lake  processes  by  major  north-flowing  riv  ?rs  and 
streams  that  have  meandered  across  this  flat  part  of  the  Coastal  Plain  following  the  last 
major  transgression.  For  that  reason,  many  lakes  in  the  area  were  mapped  in  a fluvial/ 
lacustrine  unit. 


Figure  I i.  Northern  Alaska  showing  limits  of  major  transgressions  proposed  by  previ- 
ous workers: 

A1.  Late  Sangamon  or  mid-Wisconsin  (O'Sullivan  1966) 

A.  Wisconsin  (Black  1964 ) 

B.  Sangamon  (O'Sullivan  1966) 

C.  Mid-Pleistocene  (O'Sullivan  1966) 

D.  Late  Pliocene  (O'Sullivan  1966) 

Surfaces  less  than  8 rn  and  8 to  29  m in  elevation  are  also  shown. 
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A new  delineation  of  this  transgression  is  proposed,  based  largely  on  the  inner  limit  of  a 
< 8 m surface  which  some  previous  evidence  suggests  may  have  been  the  mid-Wisconsin, 
or  prior,  sea  level  high  (Hopkins  1973).  The  use  of  elevation  data  for  the  positioning  of  this 
strand  line  assumes  that  there  has  been  no  subsequent  uplift  of  the  surface  since  the  time 
of  this  transgression;  however,  this  assumption  is  probably  not  true  for  the  entire  Coastal 
Plain.  A comparison  of  this  limit  with  previous  limits  is  in  general  agreement  except  in  the 
more  westerly  region  of  the  Coastal  Plain  where  the  limits  suggested  by  Black  (1964)  and 
O'Sullivan  (1961)  extend  farther  inland  and  westward.  The  8-m  limit  appears  to  be  more 
comparable  to  units  mapped  by  lake  morphology,  although  the  L,  class  (as  mapped)  does 
extend  farther  to  the  west  than  the  boundary  suggested  by  the  8-m  elevation.  This  evi- 
dence also  suggests  greater  uplift  in  this  area. 
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The  29-m  elevation  also  marks  a topographic  break  on  the  Coastal  Plain.  In  a restricted 
western  part  of  the  region  this  break  may  correspond  to  a pre-lllinoian  strand  line.  If  this 
line  is  carried  across  the  Coastal  Plain,  whatever  its  age  significance,  it  agrees  with  the 
transition  between  mapped  lake  units  L;l  and  L4  in  the  central  part  of  the  area.  In  the 
western  part  of  the  Coastal  Plain,  the  correlation  of  this  limit  with  lake  unit  boundaries  is 
not  as  good,  which  might  be  anticipated  because  of  the  range  of  geological  and  topo- 
graphical conditions  encountered.  Therefore,  the  best  correlation  between  mapped  lake 
units  and  strand  lines  appears  to  be  with  the  most  recent  transgression  apparent  on  the 
existing  land  surface,  as  indicated  by  the  8-m  contour.  The  agreement  would  appear  to 
improve  by  increasing  the  elevation  by  3 m in  the  western  section,  although  there  is  little 
basis  for  increasing  the  level  by  that  amount. 


CONCLUSIONS 

The  ERTS-I  imagery  provides  a regional  view  for  classifying  all  of  the  lakes  on  the 
Alaskan  Arctic  Coastal  Plain  based  on  their  size,  shape,  and  orientation.  This  classifica- 
tion resulted  in  a grouping  of  the  lakes  into  six  lake  units  and  three  nonlake  units.  A means 
of  estimating  the  depths  of  many  lakes  in  this  region  appears  possible  based  on  the 
multispectral  sequential  ERTS-I  imagery.  This  technique  relies  on  sequential  imagery  and 
observations  of  the  persistence  of  ice  cover  during  the  thaw  season.  The  lakes  that  retain 
an  ice  cover  longest  into  the  summer  are  the  deepest.  Three  lake  classes  were  developed 
based  on  this  approach;  < 1 m,  1-2  m,  and  > 2 m.  The  intermediate  depth  class  depended 
on  the  contrast  in  the  imagery  between  MSS  bands  5 and  7.  The  band  5 imagery,  with  its 
greater  potential  for  water  penetration,  made  it  possible  to  detect  lakes  that  were  frozen  to 
the  bottom  and  whose  ice  cover  was  flooded,  only  making  them  appear  ice-free. 

The  potential  maximum  depth  of  thaw  settlement,  based  on  the  volume  of  ground  ice 
and  thaw  settlement  in  the  Barrow  area,  indicated  that  conditions  on  this  part  of  the 
Coastal  Plain  favored  formation  of  lakes  less  than  2 m in  depth  and  that  lakes  deeper  than 
this  should  be  limited  in  number. 

The  relationship  between  regional  slope  and  lake  size  indicates  the  extremely  low  re- 
gional slope  required  for  development  of  large  thaw  lakes.  Most  lakes  occur  at  regional 
relief  angles  of  less  than  1°. 

The  percentage  of  land  covered  by  lakes  is  not  as  large  as  one  might  judge  from  maps, 
photos,  and  aerial  views.  The  contoured  water-land  percentages  help  to  illustrate  the 
regional  controls  of  lake  distribution.  In  addition  to  relief  one  of  the  more  pronounced 
controls  appears  to  be  the  location  of  major  drainage  networks,  which  tends  to  modify  and 
control  lake  development  and  distribution. 

The  correlation  of  the  various  mapped  lake  classes  with  land  surfaces  of  various  ages  is 
most  obvious  between  the  youngest  seaward  section  of  the  area  and  the  largest  lake  class. 
Additional  correlation  with  the  older  surfaces  is  apparent  in  some  areas  and  these  can  be 
correlated  with  both  regional  and  local  relief  patterns  unique  to  the  various  surfaces. 
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